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01 \ Abstract: 

Let S be a finite subset of a field. For multivariate polynomials the generalized 
£NJ \ Schwartz-Zippel bound [2] , [4] estimates the number of zeros over S x • • ■ x S 

counted with multiplicity. It does this in terms of the total degree, the number 
of variables and \S\. In the present work we take into account what is the 
leading monomial. This allows us to consider more general point ensembles 
and most importantly it allows us to produce much more detailed information 
' about the number of zeros of multiplicity r than can be deduced from the 

generalized Schwartz-Zippel bound. We present both upper and lower bounds. 
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^ ■ 1 Introduction 

In this paper we consider multivariate polynomials over an arbitrary field F. 
Our studies focus on the zeros of given prescribed multiplicity, a concept to 
be defined more formally below. The definition of multiplicity that we will 
use relies on the Hasse derivative. This derivative coincides with the usual 
analytic derivative in the case of polynomials over the reals. Before recalling the 
definition of the Hasse derivative let us fix some notation. Assume we are given 
a vector of variables X = (X±, . . . , X m ) and a vector k = (k±, . . . , k m ) € Nq" 
then we will write X k = X^ 1 • • • X^ 1 . We will always assume that X and Z 
are vectors of m variables. 

Definition 1. Given F(X) E F[X] and k € the k'th Hasse derivative of 
F, denoted by F^ k \X) is the coefficient of Z k in F(X + Z). In other words 



F{X + Z) = ^F^\X)Z k 



The concept of multiplicity for univariate polynomials is generalized to mul- 
tivariate polynomials in the following way. 
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Definition 2. For F{X) € FLY]\{0} and a 6 F' m we define the multiplicity 
of F at a denoted by mult(F, a) as follows. Let M be an integer such that for 
every k = (k±, . . . , k m ) <G with ki + ■ ■ ■ + k m < M, F^(a) = holds, but 

for some k = (k±, . . . , k m ) £ N™ with k\ H h k m = M, F$\a) + holds, 

then mult(F,a) = M. If F = then we define mult(F,a) = oo. 

It is of evident interest to investigate for multivariate polynomials F and a 
finite ensemble of points the following questions: 

Ql How many zeros can F have in total when counted with multiplicity? 
Q2 How many zeros of a given prescribed multiplicity can F have? 

Clearly, assuming finite ensembles of points is not a restriction when F is a 
finite field F q . We note that the above questions have important implications 
in a number of applications, see [4] and [11]. What we would like to have for 
certain natural ensembles of points is bounds on the number of points in terms 
of the total degree of F or even better in terms of lm(i^). Here, lm(i^) denotes 
the leading monomial of F with respect to some fixed monomial ordering. 
The related problem of bounding the number of zeros (counted without multi- 
plicity) has been dealt with using two completely different approaches. On the 
one hand a tight bound in terms of the leading monomial has be derived using 
the footprint bound from Grobner basis theory (see [3] and [6]). On the other 
hand a tight bound in term of the total degree, known as the Schwartz-Zippel 
bound, was derived using only very simple combinatorial arguments [12], [13]. 
To answer partly question Q2 in terms of the total degree Pellikaan and Wu 
in [11] followed the footprint bound approach. Later a generalized Schwartz- 
Zippel bound that deals with question Ql in terms of the total degree was 
suggested by Augot, El-Khamy, McEliece, Parvaresh, Stepanov, Vardy in [1] 
for the case of two variables, and by Augot, Stepanov in [2] for arbitrarily 
many variables. The bound was proven to be correct in a recent paper by Dvir, 
Kopparty, Saraf and Sudan [4]. The generalized Schwartz-Zippel bound goes 
as follows. 

Theorem 3. Let F{X) E FLY] be a non-zero polynomial of total degree d. 
Then for any finite set S C F 

mult(F,a) < d\S\ m ~ l . 

a£S n 

As a corollary we get an immediate partial answer to question Q2 in terms 
of the total degree of F. 

Corollary 4. Let F(X) £ ~F[X] be a non-zero polynomial of total degree d and 
let S C F be finite. The number of zeros of F of multiplicity at least r from S m 
is at most 

r 

In the present paper we take the Schwartz-Zippel approach. We use the 
methods from [4], but rather than taking into account only information about 
the total degree and allowing only point ensembles S n we 
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• use information about the leading monomial with respect to a lexico- 
graphic ordering. 

• consider the more general point ensembles Si x • • • x S m , the sets Si all 
being finite. 

In Section 2 we easily translate Theorem 3 into this setting and derive an imme- 
diate translation of Corollary 4. As will be shown in Section 6, Theorem 3 and 
its translation are tight for all products of univariate linear terms. A similar 
result by no means holds for Corollary 4 and its translation. Actually, a refine- 
ment of the methods from [4] yields for dramatic improvements to Corollary 4 
and its translation. In its most general form in Section 3 we state an algorithm 
to upper bound the number of zeros of multiplicity at least r. Using this algo- 
rithm we then derive in Section 4 closed formulas in the case where the number 
of variables is two and the multiplicity is arbitrary. Section 5 further presents 
a simple closed formula for the case of arbitrary many variables where, how- 
ever, the powers %i, . . . , i m in the leading monomial lm(.F) = X^ 1 • • • X % ™ are all 
small. In Section 6 we consider the case when the polynomial is a product of 
univariate linear terms. Such polynomials are easy to analyze and by doing this 
we get in appendix A an algorithm to produce lower bounds on the maximal 
attainable number of zeros of multiplicity at least r. Section 7 describes vari- 
ous conditions under which our upper bound equals our lower bound. Having 
improved on the results in [4, Section 2] we conclude the paper by showing in 
Appendix B that Corollary 4 is stronger than the corresponding result given 
by Pellikaan and Wu in [11]. From this we can conclude that the results found 
in the present paper are the strongest known. The present paper comes with a 
webpage [8] where a large number of experimental results are presented. 

2 Using information about the leading monomial 

In the following we modify the method from [4, Section 2]. One could choose to 
prove the results of the present section using the original method, however, the 
modification will be needed in the section to follow. For simplicity we stick to 
the modified method in both sections. Throughout the paper Si, ... , S m C F 
are finite sets and we write si = \Si\, . . . , s m = \S m \. In the following the mono- 
mial ordering -< on the set of monomials in variables X\ , . . . , X m will always 
be the lexicographic ordering with X m -<•••-< X\ . 

We start our investigations by recalling two results from [4, Section 2]. The 
first corresponds to [4, Lemma 5]. 

Lemma 5. Consider F(X) G F[X] and a G F m . For any k = (ki, . . . , k m ) G 
N™ we have 

mult(F^ k \a) > mult(F,a) - (h -\ h k m ). 

The next result that we recall corresponds to the last part of [4, Proposition 

6]. 
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Proposition 6. Given F(X±, . . . , X m ) € F[Xi, . . . , X m ] and 

Q(Y 1 ,...,Y l ) = (Q 1 (Y), Q m (Y)) eF[y lr ., Y t ] m 
let F o Q be the polynomial F(Q\(Y), . . . , Q m {Y)). For any a € we have 

mult(F o Q,a) > mult(F, Q(a)). 
We get the following Corollary, which is closely related to [4, Corollary 7]. 

Corollary 7. Let F(X l , X m ) € F[Xi, X m ] and h, . . . , 6 m _i, c 6 F m 
be given. Write F*(Ti, . . . , T m ,_i) = F(T\b\ + • • • + T m _i& m _i + c). For any 
(ii, . . . ,t m -i) € F m_1 we have 

mult(F*(Ti, . . . ,T m _i), (ti, . . . ,t m -i)) 

> mult(F(X 1 ,. . .,X m ),t x bi H h i m _i& m _i + c). 

We now write 

F(X 1 ,. . . , X m ) = Xf • • • X^Si Fj u ... jm _ 1 (X m ). 

Let X^ 1 • • • X^ 1 be the leading monomial of F with respect to -< . Then due to 
the definition of -<, Fi lt ...,i m _ 1 (X m ) is a (univariate) polynomial of degree i m . 
For a m £ F define 

r(a m ) = mult(F ilr .. iim _ 1 (X m ),a m ). 

Clearly, 

^ r(o m ) < i m . (1) 

We have 

^•-°^(X 1 ,...,X ro )= ^ Xf ■■■Xt-- 1 1 F^±_ 1 (X m ) 

jlvjm-l 

and due to the definition of -< and to the definition of r(a m ) we have 

l m ^( F (0,..,0,r(a m )) (Xij )Xm _ 1 ,a m )) = XI 1 • ■ • jfc 1 . (2) 

Applying first Lemma 5 with k = (0, . . . , 0, r(a m )) and afterwards Corollary 7 
with bi = (1, 0, . . . , 0), . . . , 6 m _i = (0, . . . , 0, 1, 0), c = (0, . . . , 0, a m ) and ii = 
ai, . . . ,t m -i = a m -i we get the following result which is closely related to a 
result in [4, Proof of Lemma 8]: 

mult(F(Xi, . . . ,X m ), (oi, . . . ,o m )) 

< (0 + • • • + + r{a m )) + mult( J F(°'-'°' r ^))(X 1 , . . . ,X m ), (oi, . . . , a m )) 

< r{a m ) + mult(F(°""'°' r -^))(X 1 , . . .,X m _!,a m ), (a 1; . . . ,a m _i)). (3) 

We are now in the position that we can prove the main result of this section. 
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Theorem 8. Let F(X) G ~F[X] be a non-zero polynomial and let lm(F) = 
Xj 1 • • • X^ 1 be its leading monomial with respect to a lexicographic ordering. 
Then for any finite sets S± , . . . , S m C F 

^ mult(F, a) <i\S2--- s m + siz 2 s 3 ■ • • s m H h si ■ ■ ■ s m -\i m . 

SeSix— xS m 

Proof. We prove the theorem for the monomial ordering -<. Dealing with gen- 
eral lexicographic orderings is simply a question of relabeling the variables. 
Clearly the theorem holds for m = 1. For m > 1 we consider (3). Assuming 
the theorem holds when the number of variables is smaller than m we get by 
applying (1) and (2) the following estimate 



mult (F, a) 



aeSi x— x5 m 

< ~i m Sl ■ ■ ■ S m -i + S m (iiS 2 ■ ■ ■ Sm-l H h *m-l*l • • • S m -2) 

= HS2 ■ ■ ■ S m + i2SlS 3 ■ ■ ■ S m -\ i m S\ ■ ■ ■ S m -l 

as required. □ 

We have the following immediate generalization of Corollary 4. 

Corollary 9. Let F(X) G F[X] be a non-zero polynomial and let lm(F) = 
X % y ■ ■ ■ X 1 ™ be its leading monomial with respect to a lexicographic ordering. 
Assume Si, ... , S m ^ F are finite sets. Then over Si X • • • X S m the number of 
zeros of multiplicity at least r is less than or equal to 

(iis 2 ■■■ s m + sii 2 s 3 ■■■Sm-l h si ■ ■ ■ s m -iim)/r. 

3 Improvements to Corollary 9 

In this section we shall see that a further analysis allows for dramatic 
improvements to Corollary 9. Let Xj 1 • • • XJ^ be the leading monomial of 
F{X) G F[X] with respect to -<. Recall from (3) the bound 

mult(F(Xi, . . . ,X m ), (ai, . . . ,a m )) 

< r(a m ) + mult(F( ""' ' r ( a ™))(X 1 , . . .,X m -i,a m ), (a ls . . . ,a m -i)). (4) 

Here, r(a m ) are numbers that when summed over all possible a m G S m give 
at most i m and the leading monomial of i^(°>---' > r ( am ))(Xi, . . . , X m -i, a m ) with 
respect to -< is X n ■ ■ ■ X l ™Z\ ■ Our analysis suggests the following recursive 
definition of a function to bound the number of zeros of multiplicity r. 

Definition 10. Let r G N, ii, . . . , i m G No- Define 

D(it,r,si) = min{|_— J,si} 
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and for m > 2 



D(h,... ,i m ,r,si 



max 

(ui ,. . .,u r )eA(i m ,r,s m ) 



s m -u± - ■■ 



u r )D(ii, . . . , 

^ra— 1 ; fi Si j • • • j ^m- 1 J 



+ uiD(ii, . . 



*m-i,r - 1, Si 



• • • ,s m _i) H 



+ n r _iD(ii 



• ; im—li 1 j Si , . . . , S m — 1 ) -|- U r Si ■ • • S m —1 



where 

{(ui, . . . , u r ) G Nq j u\ + • • • + u r < s m and U\ + 2u 2 + • • • + ru r < i m }. 

Theorem 11. For a polynomial F(X) G F[X] let XJ 1 ■ ■ ■ X 1 ™ be its leading 
monomial with respect to < ( this is the lexicographic ordering with X m -<•••-< 
X\). Then F has at most D(i\, . . . ,i m ,r, s±, . . . , s m ) zeros of multiplicity at 
least r in Si x • • • x S m . The corresponding recursive algorithm produces a 
number that is at most equal to the number found in Corollary 9 and is at most 
equal to Si • • • s m . 

Proof. The proof of the first part of the proposition is an induction proof. The 
result clearly holds for m = 1. Given m > 1 assume it holds for m — 1. For d = 
1, . . . , ? — 1 let Ud be the number of a m 's with r(a m ) = d and let u r be the number 
of a m 's with r{a m ) > r. The number of a m 's with r(a m ) = is s m — u\ — ■ ■ - — u r . 
The boundary conditions that u\ + ■ ■ ■ A-u r < s m and u\ + 2u 2 + • • • + ru r < i m 
are obvious. For every a m with r(a m ) = d, d = 0, . . . , r — 1 for [a\, . . . , a m ) to 
be a zero of multiplicity at least r the last expression in (4) must be at least 
r — d. For a m with r(a m ) > r all choices of a±, . . . , a m _i are legal. This proves 
the first part of the proposition. As both Corollary 9 and the above proof rely 
on (4), Theorem 11 cannot produce a number greater than what is found in 
Corollary 9. The condition u± + ■ — h u r < s m and the definition of D(i±,r, si) 
imply the last result. □ 

The next remark shows that we need only apply the algorithm to a restricted 
set of exponents (ii, . . . , i m ). 

Remark 12. Given (ii, . . . , i m , r, s±, . . . ,s m ) with H 1- [im/s m ] > r 

then there exist polynomials with the leading monomial being X[ x ■ ■ ■ X^ such 
that all points in S% X • • • X S m are zeros of multiplicity at least r. Hence, we 
need only apply the algorithm to cases that do not satisfy the above inequality. 
In Section 6, Example 31, we will explain this fact in more detail. 

In a series of experiments we found that the above algorithm produces num- 
bers that are often much lower than the minimum of the corresponding result 
from Corollary 9 and si • • • s m . In the webpage [8] we list all results of our 
experiments. Here, we only mention a few. 
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Table 1: D(ii, Z2, 3, 5, 5) 





A 

u 


i 
i 


Z 


q 


i 

4 


r. 















5 


5 


5 


1 








1 


5 


(3 


(3 


2 





1 


2 


7 


8 


9 


3 


5 


5 


5 


9 


9 


10 


4 


5 


5 


6 


9 


11 


13 


5 


5 


6 


7 


11 


12 


14 


6 


10 


10 


10 


13 


14 


17 


»2 7 


10 


10 


11 


13 


15 


18 


8 


10 


11 


12 


15 


17 


21 


9 


15 


15 


15 


17 


18 


22 


10 


15 


15 


16 


17 


20 




11 


15 


16 


17 


19 


21 




12 


20 


20 


20 


21 


22 




13 


20 


20 


21 


21 


23 




14 


20 


21 


22 


23 


24 





6 


t—r 

1 


8 


9 


10 


11 


12 


13 


1 l 

14 


10 


10 


10 


15 


15 


15 


20 


20 


20 


11 


11 


12 


1(3 


17 


17 


21 


21 


21 


13 


13 


14 


17 


19 


19 


22 


22 


22 


14 


14 


16 


18 


21 


21 


23 


23 


23 


16 


16 


18 


19 


23 


23 


24 


24 


24 


17 


17 


20 


20 












19 


19 


21 


21 












20 


20 


22 


22 












22 


22 


23 


23 












23 


23 


24 


24 













Example 13. In this example we bound the number of zeros of multiplicity 
3 or more for polynomials in two variables. Both Si and S2 are assumed to 
be of size 5. Table 1 shows information obtained from our algorithm for the 
exponents i\,i<i not treated by Remark 12. Table 2 illustrates the improvement 
on the bound Lmin{(ii + 42)5/3, 5 2 }J . Here, the first expression comes from 
Corollary 9 and the last expression is the number of points in S\ x $2- Observe, 
that the tables are not symmetric meaning that D(ii,i2, 3, 5, 5) does not always 
equal D(i2, ii, 3, 5, 5). 
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Table 2: Improvements found in Example 13 



















k 





















1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


u 


u 


1 


Q 




u 


1 




u 


1 


Q 




u 


1 


Q 

o 


U 


i 


Q 



1 

1 


1 

1 


Q 



1 

4 


1 


o 

z 


4 


U 


o 

z 


Q 




U 


i 

l 


Q 

o 


U 


n 

z 


1 

4 


n 
Z 


Q 



1 

4 


4 


1 

1 


n 

2 


o 

z 


u 


o 

z 


o 

z 


1 


1 


1 


i 

1 


Q 




Q 
O 


Q 



u 


1 


3 


1 


n 
Z 


Q 




1 

1 


n 

z 


o 

z 


o 

z 





o 
_/ 


n 
_/ 


rj 
Z 


Z 


1 

4 


1 

1 


o 



i 

4 


o 

2 


o 
z 


o 
z 


u 


o 
Z 


o 
z 


o 

z 





o 

I 


1 


1 


1 


5 


3 


4 


4 


2 


3 


2 


1 


3 


1 


3 












6 





1 


3 


2 


2 


1 


1 


2 


2 


4 












»2 7 


1 


3 


4 


3 


3 


2 


1 


3 


3 


3 












8 


3 


4 


4 


3 


3 





1 


3 


2 


2 
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1 


3 


3 


3 


1 


2 


2 


1 


1 












10 


1 


3 


4 


4 


3 






















11 


3 


4 


4 


4 


4 






















12 





1 


3 


4 


3 






















13 


1 


3 


4 


4 


2 






















14 


3 


4 


3 


2 


1 























Example 14. In i/ws example we bound the number of zeros of multiplicity 3 
or more for polynomials in four variables. The sets S\, S2, S3 and S4 are all 
assumed to be of size 6. Table 3 shows information obtained from our algorithm 
for a small sample of values (21,22^3 = 3, 24 = 5). Table 4 illustrates the 
improvement on the bound min{(ii + ^2 + ^3 + i4)6 3 /3,6 4 }. Here, the first 
expression comes from Corollary 9 and the last expression is the number of 
points in S\ x S2 x S3 x S4. 



Table 3: D(i\, ^2, 13 = 3, £4 = 5, 3, 6, 6, 6, 6) 



>2 








1 


2 


3 


4 


5 


6 


7 





468 


486 


504 


642 


666 


720 


912 


912 


1 


486 


501 


536 


651 


705 


764 


964 


964 


2 


504 


536 


574 


700 


759 


840 


1024 


1024 


3 


642 


651 


666 


771 


816 


908 


1077 


1077 


4 


666 


684 


732 


807 


880 


984 


1140 


1140 


5 


720 


750 


816 


876 


952 


1056 


1197 


1197 


6 


912 


912 


960 


976 


1024 


1134 


1260 


1260 


7 


912 


928 


980 


1008 


1060 


1155 


1263 


1263 
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Table 4: Improvements found in Example 14 



h 



12 








1 


2 


3 


4 


5 


6 


7 





108 


162 


216 


150 


198 


216 


96 


168 


1 


162 


219 


256 


213 


231 


244 


116 


188 


2 


216 


256 


290 


236 


249 


240 


128 


200 


3 


150 


213 


270 


237 


264 


244 


147 


219 


4 


198 


252 


276 


273 


272 


240 


156 


156 


5 


216 


258 


264 


276 


272 


240 


99 


99 


6 


96 


168 


192 


248 


272 


162 


36 


36 


7 


168 


224 


244 


288 


236 


141 


33 


33 



Table 5: Maximum improvements relative to q m ; truncated 
m 2 34 

r2345234523 

2 0.25 0.25 0.25 0.25 0.25 0.375 0.375 0.375 0.312 0.375 

3 0.222 0.222 0.222 0.222 0.296 0.296 0.296 0.296 0.296 0.333 

4 0.187 0.187 0.187 0.187 0.281 0.25 0.25 0.265 0.316 0.289 
q 5 0.24 0.16 0.16 0.2 0.256 0.256 0.232 0.24 0.307 0.288 

6 0.222 0.194 0.166 0.166 0.277 0.25 0.231 0.212 0.293 0.287 

7 0.204 0.204 0.163 0.142 0.279 0.244 0.227 0.209 0.299 0.276 

8 0.234 0.203 0.171 0.140 0.275 0.25 0.214 ? 0.299 0.275 



Example 15. Let s± = ■ ■ ■ = s m = q. Our experiments listed in [8] show that 
the value D{i\, . . . , i m , r,q, . . . ,q) often improves dramatically on the previous 
known bounds. We here list the maximal attained improvement for a selection 
of fixed values of m,q,r. We do this relatively to the number of points in 
Si x • • • x S m . In other words we list in Table 5 the value 

( max {min{(ii + ■ ■ ■ i m )q m ~ l /r, q m } - D(h, . . . , i m , r, q, . . . , q)} J jq m 

\il,—,im J 

for various choices of m,q,r. The experiments also show a distinct average 
improvement. This is illustrated in Table 6 where for fixed q, r, m we list the 
mean value of 

min{(n H Hm)g m ~ 1 ,q m } ~ D(k, ■ ■ ■ ,i m ,r,q, . . . ,q) 

mm{(i 1 + ---+i m )q m " 1 ,q m } ' 1 ' 

The average is taken over the set of exponents (i\, . . . ,i m ) / where [h/q\ + 
• • • + \jm/l\ < t holds. 

Example 16. In Example 13 for any total degree d there exists a choice of 
ii, . . . , i m with i\ + • • • + i m = d such that 

D(h, ■ ■ ■ ,i m ,r,q, ... ,q) = min{dg m_1 /r,q m }. 
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Table 6: The mean value of (5); truncated 



m 2 3 4 
r2345234523 

2 0.363 0.273 0.337 0.291 0.301 0.300 0.342 0.307 0.248 0.260 

3 0.217 0.286 0.228 0.236 0.194 0.224 0.213 0.214 0.158 0.177 

4 0.191 0.197 0.232 0.195 0.158 0.169 0.180 0.172 0.125 0.135 
q 5 0.155 0.167 0.174 0.197 0.139 0.145 0.148 0.153 0.110 0.116 

6 0.148 0.160 0.156 0.154 0.128 0.132 0.132 0.131 0.100 0.105 

7 0.128 0.137 0.138 0.138 0.119 0.122 0.121 0.119 0.093 0.098 

8 0.126 0.127 0.134 0.126 0.114 0.115 0.113 ? 0.089 0.093 



However, there are cases where such a result does not hold. Going through all 
possible choices of with ^i + «2 + *3 = 12 we see that the largest obtained 

value of D{i\, i2, is, 3, 8, 8, 8) equals 224 whereas min{12 • 8 2 /3, 8 3 } = 256. 

The next two examples are of a theoretical nature. 

Example 17. Given an arbitrary monomial ordering let lm(F) = X 1 ^ ■ ■ ■ X 1 ™ 
with i\ < si,...,i m < s m . Using results from Grobner basis theory we can 
deduce that F can have no more than 

si ■ ■ ■ s m - (si - ix) ■ ■ ■ (s m - i m ) (6) 

zeros (of multiplicity 1 or more) over S\ x • • • x S m . (See [7] and [5] for the case 
of S\ = ■ ■ ■ = S m = F q .) This result is known to be sharp meaning that polyno- 
mials exist with this many zeros. It is interesting to observe that (6) follows as 
an immediate corollary to Theorem 11 in the case where the monomial ordering 
-< is the pure lexicographic ordering with X m -<•••-< X\ . In contrast (6) only 
equals the result in Corollary 9 when lm{F) is univariate; in general the two 
bounds can differ very much. In Section 5 we will see that for the case of the 
monomial ordering being -< (6) can be viewed as a special case of a more general 
result. 

Example 18. Consider that the leading monomial is univariate, i.e. lm(F) = 
X l t for some t 6 {1, . . . , m}. Theorem 11 tells us that F can have at most 

i i 

si • •• st-i[_—\s t +i ■■■s m 
r 

zeros of multiplicity r or more. In contrast Corollary 9 only gives us the bound 

i i 

[Sl ■ ■ ■ St-l — St+l ■ ■ ■ S m \- 

r 

For m > 1 the bounds are the same only when r divides it- Assume in larger 
generality that it i: ■ ■ ■ ,it v , t u < t w for u < w are the non-zero elements in 
{ii,...,i m }. Then 

D{i 1 ,...,i m ,r,s 1 ,...,s m ) = ( s d J D(i tl , . . . , i tv , r, s tl , . . . , s tv ). 
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4 The case of two variables 



In this section we derive closed formulas for the case of two variables and the 
multiplicity being arbitrary. By Remark 12 the following Proposition covers all 
non-trivial cases. 

Proposition 19. For k = 1, . . . , r — 1, D{ki i 2 , r, si, s 2 ) is upper bounded by 

(C.l) S 2 — + 



r r r — k 

if (r — k)j^s\ < i\ < (r — k)s\ and < i 2 < ks2 

(C.2) s 2 ^ + ((k + l)s 2 - k){-^T --) + (»2 " ks 2 )( Sl - ^) 
r r — k r r 

if (r - < n < (r - A;)si and /es 2 < «2 < (k + l)s 2 

r k + 1 r 
if (r — k — l)s\ < i\ < (r — A;)^-si and < z 2 < (k + l)s2- 

^1 ^ 1 

fC.^J D(h,i 2 ,r, s 1 ,s 2 ) = s 2 [—\ + i 2 (si - |_— J) 

r r 

i/ si(r — 1) < i\ < sir and < i 2 < s 2 . 
The above numbers are at most equal to min{(iis 2 + s\i 2 )/r, s\s 2 }. 

Proof. First we consider the values of k, i 2 , r, si, s 2 corresponding to one of the 
cases (C.l), (C.2), (C.3). Let k be the largest number (as in Proposition 19) 
such that k < (r — k)s±. Indeed k G {1, . . . , r — 1}. We have 

D(i 1 ,i2,r,s 1 ,S2) < 

f k . , k k, . . , k k, 

max \s 2 \-ui{ )H h u k ( ) 

(ui,...,ur)eB(i 2 ,r,s 2 ) l r r — 1 r r — k r 

+ u fc+ i(si - — ) H Vu r {si-— )\ (7) 

r r J 

where 

B(i 2 ,r,s 2 ) = . ,u r ) € Q r I < u u . . . ,u r , 

u\ + h u r < s 2 ,ui + 2u 2 + h ru r < i 2 }. 



We observe, that 



H^i--)<K- H 11 



'r — l r r — k r 

holds for I < k. Furthermore, we have the biimplication 

(r - fc)-r7«i < k O (fc + l)(^-r - -) > - -). 
r+1 r — k r r 

Therefore, if the conditions in (C.l) are satisfied then (7) takes on its maximum 
when Uk = 77 an d the remaining itj's equal 0. If the conditions in (C.2) are 
satisfied then (7) takes on its maximum at u k = (k+ l)s 2 — i 2 , Uk+i = (i 2 — ks 2 ) 
and the remaining u^s equal 0. If the conditions in (C.3) are satisfied then (7) 
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takes on its maximal value at Uk+i = and the remaining u^s equal 0. 
Finally, if si(r — 1) < i\ < sir and < 12 < S2 then D{i\, i%, r, s±, S2) is the 
maximal value of 

S2 L — J + - L-J) + ■ • ■ + «r(«l " L-J) 

over B(i2,r, 82)- The maximum is attained for ui = 22 an d all other Uj's equal 0. 
The proof of the last result follows the proof of the last part of Theorem 11. □ 

Remark 20. Experiments show (see [8]) that the numbers produced by Propo- 
sition 19 are often much smaller than min{(iiS2 + s\i2)/r, s\S2}- However, 
there are cases where they are the identical. This happens for example when 
h = si(r — 1) and r divides s\ and S2- In the proof of (C.l), (C.2), (C.3) 
we allowed Ui,...,u r to be rational numbers rather than integers. Therefore 
we cannot expect the upper bounds in Proposition 19 to equal the true value of 
D(i\, 12, r, s\, S2) in general. Our experiments show that the bounds in (C.l), 
( C.2), ( C.3) are sometimes close to D(i±,i2, r, s±, S2) but not always. Hence the 
best information is found by actually applying the algorithm from the previous 
section. 



5 When ii, . . . , i m are small 

Having already four different cases when m = 2 the situation gets rather com- 
plicated when we have more variables. Assuming however that all exponents 
i\, . . . , im in the leading monomial are small we can give a very simple formula. 
Whereas the formula is simple we must admit that the precise definition of 
i±, . . . , i m being small is a little involved. It goes as follows. 
Definition 21. Let m > 2. We say that . . . , i m ,r, s±, . . . , s m ) satisfies 
Condition A if the following hold 

(A.l) ii,...,i m <s m 

(A2) s ( Sl -k)...( Sm _ 2 _in^)< Z ( Sl _i) 

for all I = 2, . . . , r, s = 1, . . . I — 1. 
(A3) s ( Sl -^)...(s m ^-^)<r( Sl -^) 

for all s = 1, . . . , r — 1. 

Example 22. If r = 1 then (A. 2) and (A. 3) do not apply and with a reference 
to Remark 12 (A.l) is a natural requirement. 

Example 23. For m = 2 and r arbitrary condition (A. 2) does not apply and 
condition (A. 3) simplifies to 

2 2 
r s — rs 

h < — n y-si 

r z — s z 

for all s with 1 < s < r. The minimal upper bound on i\ is attained for s = 1. 
Hence, in case of two variables Condition A reads i\ < fxfSi, i2 < S2- 



■■■{Sm-2-^) 
■ ■ ■ {Sm-l ~ ^) 
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2/3 



Figure 1: The surface § {h + h + h) - \{hh + hh + hh) + fhhh = 1 
Example 24. For r = 2 conditions (A. 2), (A. 3) simplifies all together to 

( s l - 4) - - " ( s m-l - < 2 ( s l - *l) ■ ■ ■ - im-l)- 



For r = 2, m = 3 and si = S2 = S3 = q Condition A therefore reads 

^(I 1+ I 2 )- 7 -I 1 I 2 <l, / 3 <1 

where I\ = ii/g, I2 = 22/0 and F3 = 23/0. For r = 2, m = 4 and s\ = si = 
S3 = S4 = q Condition A reads 

3 7 15 

-{h + h + h) - ^(hh + /1/3 + hh) + -jhhh < 1, ^4 < 1 

where I4 = i^/q. This is illustrated in Figure 1. 

In Example 17 we discussed a well known bound on the number of zeros 
of multiplicity at least r = 1. With Example 22 in mind the last part of the 
following Proposition can be viewed as a generalization of this bound. 

Proposition 25. Assume that (ii, . . . , i m , r, s\, . . . ,s m ), m > 2 satisfies Con- 
dition A. If r > 2 then 

r r 
k< — — 7S1, . . . ,i m -i < — —Sm-i- (8) 
r + 1 r + 1 
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For general r we have 



D(i 1 , . . . , i m , r, si, . . . , s m ) < si ••• s m — (si -—)•■■ (s m - — ) (9) 

which is at most equal to min{(i 1 S2 • • • s m + ■ ■ ■ + S\ ■ ■ ■ s m -ii m ) /r, s\ ■ ■ ■ s m }. 
Proof. We start by noting that (A. 2) implies 

(si - j) ■ ■ ■ (s t -i - < l(si - —) - ■ ■ (s t -i - 

for all i = 2, . . . , m — 1, I = 2, . . . , r, s = 1, . . . I — 1. A similar thing holds 
regarding (A. 3) and if we combine this fact with the result in Example 23 we 
get (8) for r > 2. Now let . . . ,i m , r, s±, . . . , s m ) with m > 1 be such that 
Condition A holds. We give an induction proof that 

D(h, . . . ,i t , I, si, . . . , s t ) < si ■ • • s t - (si - ^) • • • (s t - ^) (w) 

for all 1 < t < m, 1 < I < r { ' 

For t = 1 the result is clear. Let 1 < t < m and assume the result holds when 
t is substituted with t — 1. According to Definition 10 we have 

D(ii,...,i t ,l,si,...,s t ) = 



max <( s t - ui ui)D(h, . . .,i t -i,l,S!, . . . ,s t -i) 

(ui,...,ui)GA(it,l,st) I, 

+ uiD(n, it-x,l - 1, si, . . . , s t _i) H 

+ ui-iD(h, . . -,it-i, l,si, . . . ,s t -i) + uisi ■ ■ ■ s t -i 

where 

A(i t ,l,s t ) = {(ui,...,ui) G No | u\ +2u 2 -\ \- lui < i t } 

follows from (8). By the above assumptions this implies that 

D(ii,...,i u l,si,...,s t ) < 



max 



s t (si ■ ■ ■ s t -i - (si - -j-) ■ ■ ■ ( s t-i ~ "V")) 



(ui,...,ui)eB(it,l,s t ) [ / / 

+ «i((ai - y) • • • ( s *-i - -y 1 ) - ( s i - r~[) ■ ■ ■ ( s *-i - frj)) 
+ ••• 

+ ((si - y) • • ■ ( s <-i - -y 1 ) - («i - y) • • • (st-i - -y^)) 
+ «i((ai - y)---(st-i - ^p)) 
where 

B (it, I, s t ) = {(ui,...,ui) eQ l \ 0<ui,...,ui and ui + 2« 2 H h /«j < it}. 
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As t < m condition (A. 2) applies and tells us that the maximal value is attained 
for u\ = ■ ■ ■ = ui-\ = and ui = ¥. This concludes the induction proof of (10). 
To show (9) we apply similar arguments to the case t = m but use condition 
(A. 3) rather than condition (A. 2). The proof of the last result in the proposition 
follows the proof of the last part of Theorem 11. □ 

Remark 26. Experiments show (see [8]) that the bound in Theorem 11 is very 
often much better than min{(iiS2 • • ■ s m +- ■ - + si • • • s m -ii m )/r, si • • • s m }, how- 
ever, they also reveal that in many cases one can get more information about 
the number of zeros by actually applying the algorithm from Section 3. 

Example 27. This is a continuation of Example 23 where we translated Con- 
dition A into bounds on i\ and i<i in the case of two variables. Applying in turn 
Proposition 25 and (C.3) in Proposition 19 with k = r — 1 we see that the two 
bounds produce the very same values when m = 2. 



6 Products of univariate linear terms 

In this section we study the situation where F(X) is a product of univariate 
linear terms. First we note that equivalently to Defintion 2 one can define the 
multiplicity of a polynomial as follows. 

Definition 28. Let F(X) G F[A]\{0} and a = (ai, ...,a m ) € F m . Consider 
the ideal 

J t = ((Ax - ai ) Pl • • • (X m - a m ) p ™ \pi + ---+p m = t) C F[X U . . .,X m }. 
We have mult(F,a) = r if F £ J r \J r+ \. If F = we have mult(F,a) = oo. 

The above definition makes it particularly simple to calculate the number 
of zeros of multiplicity at least r when F is a product of univariate linear terms. 
In the following write 

*j \«r < : \ 

for j = 1, . . . , m. 
Proposition 29. Consider 

F(X) = lHl(X u -a{? 

u=l v=l 
.,(!) J"")- 



The multiplicity of • • • , a j m ) i- n F(X) equals 



r) H h r) '. (11) 

Proof. Without loss of generality assume j\ = ■ ■ ■ = j m = 1. Clearly, the 

(1) (m) 

multiplicity is greater than or equal to r = r\ + • ■ ■ + r\ . Using Grdbner 

d) 

basis theory we now show that it is not larger. We substitute Xi = Xi — a\ 
for i = 1, . . . , m and observe that by Buchberger's S-pair criteria 

B = {XI' ■ ■ ■ X%» | n + • • • + r m = r + 1} 
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is a Grobner basis (with respect to any fixed monomial ordering). The support 
of F(X\, . . . , X m ) contains a monomial of the form X^ 1 ■ ■ ■ X^ with i\ + ■ ■ ■ + 
im = r - Therefore the remainder of F(X\, . . . , X m ) modulo B is non-zero. It is 
well known that if a polynomial is reduced modulo a Grobner basis then the 
remainder is zero if and only if it belongs to the ideal generated by the elements 
in the basis. □ 

We now show that Theorem 8 is tight. It follows of course that so is Theo- 
rem 3 (a fact that has not been stated in the literature). 

Proposition 30. Let Sx,...,S m C F be finite sets. If F(X) £ F[X] is a 
product of univariate linear factors then the number of zeros of F counted with 
multiplicity reach the generalized Schwartz- Zippel bound (Theorem 8). 

Proof. Consider the polynomial 

_ m f u _ (u) 

w=nn(^-«iT • 

U=l 11=1 

Write i u = Ylv=l r v , u = 1, . . . ,m. We have 

■S 1 Sjn 

mult(F, a) = ^(s 2 • • • s m )r, (1) + • • • + 

a£SiX---xSm * =1 t=l 

— ^1^2 " " " S m T • • ■ T Si 1 " S m —\% m 

Example 31. Let . . . , i m ,r, s±, . . . , s m ) be such that [ii/sij-h ■ •+[i m /s m \ > 
r. As mentioned in Remark 12 there exist polynomials with the leading mono- 
mial being X[ x ■ ■ ■ such that all points in S\X ■ ■ • x S m are zeros of multiplic- 
ity at least r. To see this define r\ = \ i\/ s\\, ... ,r m = [im/s m \. Multiplying 

m s u 

nn(^-4 u) r 

u=l v=l 

by an appropriate monomial we get a polynomial having the prescribed leading 
monomial (with respect to any monomial ordering). Clearly, all points in the 
ensemble are zeros of multiplicity at least r. 

Definition 32. Given . . . , i m , r, s\, . . . , s m ) consider the set of polynomials 
that are products of univariate linear terms and have X 1 ^ ■ ■ ■ X 1 ™ as leading 
monomial. By H(i%, . . . , i m ,r, s±, . . . , s m ) we denote the maximal number of 
zeros of multiplicity at least r that a polynomial from the above set can have 
over S\ x • • • x S m . 

Based on Proposition 29 it is straightforward to describe an iterative algo- 
rithm that finds H(ii, . . . , i m , r, si, . . . , s m ). For the convenience of the reader 
we include such an algorithm in Appendix A. 



-i)r 



(m) 



□ 
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Table 7: Difference between upper and lower bound in Example 33 
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In the previous sections we considered the general set of polynomials F with 
lm^(F) = X^ 1 ■ ■ ■ X 1 ™. We gave upper bounds on the maximal attainable num- 
ber of zeros of multiplicity r or more. It is clear that H(i\, . . . , i m , r, s x , . . . , s m ) 
serves as a lower bound on the maximal attainable number of zeros of multiplic- 
ity r or more. In particular H(i±, . . . , i m , r, si, . . . , s m 

1"mi r i . . . , S m ) 

holds. Experiments show (see [8]) that the two functions are sometimes quite 
close. In the next section we present various conditions under which the two 
functions attain the same value. Clearly, when this happens we know what is 
the maximal number of zeros of multiplicity at least r that any polynomial with 
leading monomial Xj 1 • • • Xffi can have over Si x • • • x S m . 

Example 33. This is a continuation of Example 13 where we studied the upper 
bound D(ii, 12, 3, 5, 5) for relevant choices o/ii,*2- Using the algorithm in Ap- 
pendix A we calculated the corresponding values of the lower bound H(ii,i2, 3, 5, 5). 
We list in Table 7 the difference D(ii,i2, 3, 5, 5) — H(ii,i2, 3, 5, 5). We see that 
for many choices o/«i,«2 the upper bound equals the lower bound. 

Example 34. This is a continuation of Example 14 where we studied D(ii,i2,i3 = 
3, ii = 5, 3, 6, 6, 6, 6, ) for a collection of values %i,%2- In Table 8 we list the dif- 
ference between these upper bounds and the lower bounds JT(ii,«2)*3 = 3, 14 = 
5, 3, 6, 6, 6, 6). 

Example 35. Let s\ = ■ ■ ■ = s m = q. Our experiments listed in [8] show that 
D(ii, . . . , i m ,r, q, . . . ,q) is often close to H(ii, . . . , i m , r,q,...,q). In Table 9 we 
list the mean value of 

D{ii, ...,i m ,r,q,...,q)- H(ii, . . . , i m , r, q, . . . , q) 
\ {D(k ,...,i m ,r,q,...,q) + H(i x , . . . , i m , r, q, . . . , q)) 
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Table 8: Difference between upper and lower bound in Example 34 
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70 
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90 
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106 
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114 
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5 
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129 
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82 


81 
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84 


6 


216 
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201 


112 


52 


54 


72 


54 


7 


186 


181 


146 


81 


40 


42 


57 


42 



The average is taken over the set of exponents with \i\/q\ + • • • + \}m/(l\ < r 
and D(ix, . . . , i m ,r, q, . . . , q) ^ 0. 

Table 9: Mean value of (12); rounded up 
m 2 3 4 

r2345234523 

2 0.044 0.066 0.08 0.088 0.048 0.085 0.106 0.120 0.041 0.081 

3 0.039 0.048 0.068 0.075 0.046 0.067 0.092 0.104 0.038 0.064 

4 0.044 0.057 0.054 0.067 0.049 0.075 0.083 0.098 0.037 0.068 
q 5 0.042 0.057 0.061 0.060 0.045 0.073 0.086 0.092 0.034 0.064 

6 0.041 0.057 0.065 0.066 0.043 0.072 0.087 0.095 0.031 0.061 

7 0.040 0.054 0.063 0.066 0.042 0.069 0.085 0.094 0.030 0.058 

8 0.038 0.053 0.062 0.067 0.040 0.067 0.083 ? 0.029 0.056 



7 Some conditions for H = D to hold 

As the polynomial ring in one variable is a unique factorization domain we 
get H(ii,r,si) = D(ii,r, s m ) for all choices of ii,r, s%. Experiments suggest 
(see [8]) that for two variables we have a similar equality for certain systematic 
choices of ii,i2- For other choices of i\,i 2 the picture is more blurred. The 
results of our experiments further suggest that it might not be an easy task to 
say much about which values of . . . -i m , r, s±, . . . , s m ) causes equality when 
m > 3. We summarize our findings below. 

Proposition 36. For ^j-si < ii < s\, (r — l)s2 < «2 < r S2 we have 

H(h,i 2 ,r, s 1 ,s 2 ) = D(i 1 ,i 2 ,r, s 1 ,s 2 ) = rs 2 i x + i 2 si - i\i 2 - (r - l)sis 2 - 
Proof. The value of D is upper bounded by (C.2) in Proposition 19. The value 
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of H is lower bounded by studying the zeros of 

(x 2 - a ®y ■ ■■(x 2 - a®y(x 2 - ogkr 1 • • • 

(X 2 -a®y-\X 1 - a W).--(X 1 -aW) 

where w = i 2 — (r — l)s 2 . d 

We leave the proofs of the following two results for the reader. 

Proposition 37. Assume r < s±. If < i\ < r and < i 2 < rs 2 holds then 

H(ii,i 2 ,r,s 1 ,s 2 ) = D(h,i 2 ,r, si, s 2 ) = [i 2 /r\s 2 + 5 

where 5 = i\ — (r — w) + 1 if r — w < %\ and 5 = otherwise. 

Proposition 38. If H(h, . . . ,i m ,r,si, . . . , s m ) = D(h, . . . ,i m ,r, si, . . . , s m ) 
then 

H(h,.. . ,i t ,0,k+i, ■ ■ ■ ,i m ,r,si, . . . ,s t ,s',s t +i, ...,s m ) 
= D(h, . . . ,i t ,0,i t+1 , . . . ,i m ,r,si, . . . ,s t ,s',s t+ i, ...,s m ) 
— s D(i\ j • • • j im > si, . . . , s m ) . 

Proposition 39. Assume (i 

I; • • • j imi fj si, . . . , s m ) satisfies Condition A (Def- 
inition 21 ) and that r divides i\, . . . ,i m . Then 

H(i\ , . . . , i m , 7", s± , . . . , s m ) — D{i\ , . . . , imi si, . . . , s m ) 

= si ■ ■ ■ s m - [Si J • • • (s m J. 

r r 

Proof. Consider 

and apply Proposition 25. □ 
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A An algorithm to calculate H 

We here give the details of the algorithm mentioned in Section 6. 
Definition 40. Consider vectors 

^ = rf\ • • .,vW), ...M m) = (v^,- ■ -A m) ) € N5. 
Let si, . . . , s m 6 N. Define for k = 1, . . . , r 
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and for k < r and m > 2 

H(v^,...M m) ,k, Sl ,...,s m ) = 

[s m - (v { r ] + ■■■ + vt ] )] ■ ■ , v ( - m ~ 1 \k, Sl ,..., Sm _x) 

+ v[ m) H(iP-\. . . , v^ m ~ l \k - 1, ai , . . . , Sm _!) 
+ • • • + v ( ™\h(^ v( m ~V , 1, s u . . . , Sm _i) 
+ J ff(^ m \fc,s m )si---s m _i. 

Proposition 41. 

H(i\, . . . , z m , r, si , . . . , s m ,) — 

max|F(t; (1) , . . . ,tT (m) ,r,si, . . . ,s 

vP + 2vf + ■■■ + rv® =i t , for t = l,...,m 



vf ] + ■■■ + 4* } < at, 



B Comparison of Theorem 3 to a bound by Pellikaan 
and Wu 

As mentioned in the introduction for S\ = ■ ■ ■ = S m = F q there is an alterna- 
tive to Dvir et al.'s method from [4], namely the method by Pellikaan and Wu 
in [10] and [11]. We conclude the paper by showing that this other approach 
is never better than Corollary 4. Thus the results in the present paper are the 
best known results. 

In [10] Pellikaan and Wu presented two algorithms for decoding general- 
ized Reed-Muller codes. The first algorithm is based on the fact that gener- 
alized Reed-Muller codes can be viewed as subfield subcodes of Reed-Solomon 
codes whereas the second algorithm is a straightforward generalization of the 
Guruswami-Sudan decoding algorithm in [9]. The analysis of the second algo- 
rithm in [10] relies on a generalization of the footprint bound from Grobner 
basis theory. As the first algorithm outperforms the second, the details of the 
analysis of the second are not included in the journal paper [10] but can be 
found in [11]. To state the generalization of the footprint bound we will need 
the following two lemmas corresponding to [11, Lemma 2.4] respectively [11, 
Lemma 2.5]. 

Lemma 42. Given a polynomial F(X) £ F 9 [X] consider the ideal 

I(q, r, m, F) = (F) + ((Xf - ■ ■ ■ (X^ - X m )^ \ ei + • ■ ■ e m = r). 

If t is the number of points in Fj? where F has at least multiplicity r, then 

(m ~\~ r — 1 \ 
r-i r 
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Lemma 43. Let d be the total degree of F(X) £ F q [X] and define w = \_d/q\. 
If d < qr then an upper bound for the dimension of 



F q [X]/I(q,r,m,f) 

is given by 

m + r — l\ m , , , x fm + r — w — 2\ m _ l (m + r — w — l\ m 



+ {d-qw){ 7 K" ~ )q n - 

m J \ m — 1 J \ m J 

Combining the two lemmas above we get the following result which is used 
in [11] without being stated explicitly. 

Proposition 44. Let the notation be as in the above lemmas and assume d < 
qr. The number of points in where F has at least multiplicity r is at most 
equal to 



m+r— 1\ m i / j „„ \ /m+r-«i-2\ m-1 fm+r— w— 1\ m 



ri(g,r,m,d) 



in 



\ r-1 I 



Augot and Stepanov [2] gave another interpretation of Pellikaan and Wu's 
second decoding algorithm in [10] by using Theorem 3 instead of Proposition 44. 
Doing this they were able to correct much more errors which indicates that the 
generalized Schwartz-Zippel bound is stronger than Proposition 44. We here 
provide a direct proof of this fact. 

Proposition 45. Let F2(q,r,m, d) =dq m ~ l /r, then 

Txiq, r, m, d) > T 2 (q, r, m, d) 
holds for all d G [0, rq — 1]. 

Proof. We consider Ti and T 2 as functions in d on the interval [0,rg]. Our 
first observation is that Ti is a continuously piecewise linear function, each 
piece corresponding to a particular value of w. The corresponding w slopes 
constitute a decreasing sequence. Combining this observation with the fact 
that T 2 is linear in d and with the fact that 

T\(q, r, m, 0) = T 2 {q, r, m, 0) and T\(q, r, m, rq) = T 2 (q, r, m, rq) 

proves the proposition. □ 
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